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Abstract 

The one-loop QCD heavy quark potential is computed to order v 2 in the 
color singlet and octet channels. Several errors in the previous literature 
are corrected. To be consistent with the velocity power counting, the full 
dependence on |p + p'|/|p' — p| is kept. The matching conditions for the 
NRQCD one-loop potential are computed by comparing the QCD calculation 
with that in the effective theory. The graphs in the effective theory are also 
compared to terms from the hard, soft, potential, and ultrasoft regimes in 
the threshold expansion. The issue of off-shell versus on-shell matching and 
gauge dependence is discussed in detail for the l/(m|k|) term in the potential. 
Matching on-shell gives a f/(m|k|) potential that is gauge independent and 
does not vanish for QED. 
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I. INTRODUCTION 



For processes involving a non-relativistic heavy quark and antiquark, it is useful to 
combine the a s expansion of QCD with an expansion in powers of the relative velocity v. 
The scattering of a heavy quark and antiquark, Q(p) + Q(— p) — > Q(p') + Q(—p'), can be 
described using a potential V, which has an expansion in powers of the velocity v and a s , 

oo oo 

^(p, p ') = v(n \ v(n) = Y. v(n ' j) ' 

n=-2 j=l 

where V (n) ~ v n , V {n ' j) ~ v n a{ . (1) 

An important complication of non-relativistic scattering is the presence of two low-energy 
scales, mv and mv 2 , which correspond to the momentum and energy of the heavy quark, re- 
spectively. In this paper, we will assume that mv 2 3> Aq C d, and ignore any non-perturbative 
effects. 

The first term in Eq. (JJ), V^~ 2,1 \ is the static Coulomb potential at tree- level, 

= C^^- 2 , (2) 
|P-P'I 

where the color factor C depends on the relative color state of the incident quark and 
antiquark. For j > 1, V^~ 2 '^ are perturbative corrections to the Coulomb potential which 
are known to two loops 

The static potential for a color singlet QQ pair is often defined in terms of the expectation 
value of a rectangular Wilson loop of width R and length T: 

V { - 2) (R) = - \im — In ^Tr P exp (ig j dx^T 11 ) \ . (3) 

The Feynman diagrams corresponding to this expectation value build up the exponential 
of the static potential. As a result, in computing the static potential V^~ 2 ^ at order j, 
iterations of lower order terms in the potential, V^~ 2 ' n \ n < j, are subtracted. At three loops 
Appelquist, Dine and Muzinich [[J have shown that infrared divergences are encountered 
in V"* -2 ' 4 -* (R) which are not canceled by simply subtracting iterations of the lower order 
potentials. Thus, using the definition in Eq. (0) at higher orders in perturbation theory, or 
generalizing this approach to subleading terms in the v expansion, becomes cumbersome 
as the set of subtractions become more complicated, and perturbative subtractions are 
insufficient to render the potential in Eq. well-defined. 

A convenient framework for investigating the expansion in Eq. (JT|) is the effective field 
theory for non-relativistic QCD (NRQCD), formulated with a power counting in v PHT6 . 
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In the effective field theory, it is more convenient to define the potential as the Wilson 
coefficient of a four-quark operator JO] 

£ P = ~ V oI3Xt(P, p') \F [VVaC^) ^ppi*) X-p'{( X ) X-pr(^)] • (4) 
P,P' 

In Eq. (PD|) , the fields ip and \ annihilate a quark and an antiquark, respectively. The 
fields are labelled by a momentum p, and a greek index for their color and spin. The 
operator in Eq. (Q) is local on the scale x ~ 1/mv 2 , but non-local on the scale p ~ mv. 
The potential V is computed as a matching coefficient at the scale p. = m between QCD 
and an effective theory for non-relativistic QCD valid below the scale m. The effective 
theory is constructed to have the same infrared structure as perturbative QCD for the two 
heavy quark system. Therefore, defining the potential as a matching coefficient provides an 
infrared safe definition. For instance, Ref. |17| shows how the three-loop matching potential 
V is infrared safe, despite the divergence in the QCD potential of Appelquist, Dine and 
Muzinich @. 

Although several different formulations of the effective theory for non-relativistic QCD 
are currently in use [p|- p^6| , p^ -p0| , certain universal features have emerged. The on-shell de- 
grees of freedom in the effective theory include quarks with energy E ~ mv 2 and momentum 
p ~ mv, soft gluons with E ~ p ~ mv, and ultrasoft gluons with E ~ p ~ mv 2 . The soft 
and ultrasoft modes are distinct, for instance a consistent power counting in v demands that 
the ultrasoft gluon interactions are multipole expanded while soft gluon interactions 
are not. The soft gluons are essential to correctly reproduce the beta function in the ef- 
fective theory fl4] , and run between the scales m and mv. Other massless on-shell fields, 
such as light quarks, will have ultrasoft and soft components too. There are also important 
off-shell field components, such as the exchange of gluons with E ~ mv 2 and p ~ mv that 
build up the potential. Soft heavy quarks with E ~ p ~ mv are also off-shell relative to 
the heavy quark states of interest. These off-shell field components can be integrated out 
of the Lagrangian in the effective theory. Doing this leaves a Lagrangian that is non-local 
at the scale mv, but local at mv 2 . This procedure, which treats the potential components 
as four quark operators, was first seriously investigated in Ref. [PH , and the resulting ef- 



fective theory is referred to as pNRQCD. In Ref. JTTJ it was proposed that the matching 
onto effective theories should take place in two stages: at p, = m one matches QCD onto 
NRQCD as originally defined in Ref. ||, and then matches NRQCD onto pNRQCD at the 
scale /i = mv. 

The matching of four-quark operators at m was considered in Ref. JT2J, following the 
proposal in Ref. |§ that the matching procedure should be similar to that in HQET. How- 
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ever, in general this procedure seems to be slightly problematic [llj. First note that it is 
necessary to include the kinetic term for the potential quarks immediately below the scale m. 
For instance, this is necessary to reproduce the threshold value of the two-loop anomalous 
dimension for the heavy quark production current [21] in the effective theory [TjJ. With 
the kinetic energy term included at leading order, the consistency of the v power counting 
forces us to have both ultrasoft and soft modes right below m (along with the multipole 
expansion etc.). If we note that the kinetic energy term is necessary to correctly reproduce 
the Coulombic infrared divergences in the effective theory, then it might seem obvious that 
it must be included in order to have the correct infrared structure. At one-loop, matching 
exactly at threshold in dimensional regularization seemed to provide a method of avoiding 
this fll2f , but this procedure fails at two- loops. 

Since we must include the kinetic energy term for all fi < m it seems quite natural 
to immediately consider matching at \i — m onto an effective theory where the off-shell 
potential gluons and soft quarks have been integrated out. Such a formulation was proposed 
in Ref. |16| and will be considered in this paper. It will be referred to as vNRQCD. In 
this theory it is more appropriate to consider the running from m to the scales mv and 



mv in a single step. This is accomplished by using the velocity renormalization group |T6 



Once the vNRQCD effective theory has been run down to the low scale the soft degrees of 
freedom have served their purpose and may be integrated out. In this final effective theory 
no additional running needs to be considered. 

In this paper we compute the one-loop matching for the QQ and QQ potentials between 
QCD and vNRQCD up to corrections of order v 2 . Ours results are compared to the match- 



ing calculation of Pineda and Soto |12| for the four-quark operators and to the threshold 
expansion [pL3 1 . The main difference between the the non-relativistic theories in Refs. [16| 



and |TT],|12|,|T9[ is the way in which large logarithms of v are summed in the effective theory. 
There is also some difference in the matching corrections; some contributions in vNRQCD 
that arise at the scale \i = m are instead computed at /i = mv in pNRQCD. We will discuss 
the differences between the two approaches in the text of the paper. 

The full QCD calculation of the QQ scattering to order a 2 v° in the non-relativistic 
expansion has been done before p2| , p3f . In calculating the potentials, Ref. |22| performs an 
additional expansion, assuming 

(p'+P^ 2 



(P' -P) 2 



« 1 . (5) 



In the usual v power counting, p and p' are both of order v, so the ratio in Eq. (^) is 
of order unity, and cannot be treated as small. For this reason, in section II we redo the 
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QCD calculation keeping the full dependence on this ratio for both the color singlet and octet 
channels. For the order v° spin-independent potential, terms proportional to (p + p') 2 which 
were dropped in previous calculations p2|j23|| are necessary to match infrared divergences 
between the full and effective theories. 

In section III we discuss the order v° potential of the form 

(P ' 2 - p2 » 2 (6, 



171 (P' ~~ P ) 

For free quark states this potential vanishes on-shell by energy conservation, p' 2 = p 2 , and 
need not be included in the potential if one uses on-shell matching P^j . The potential in 
Eq. (P) gives a non-zero contribution in loops or in matrix elements with external Coulomb 
states, and is often included in the quark potential. For instance, the usual Breit Hamiltonian 
includes a potential of the form in Eq. (||). A potential that is only non-zero off-shell can 
be gauge dependent, and including the potential of Eq. (|) induces gauge dependence in 
the coefficients of on-shell potentials. Using an off-shell potential gives correct results if all 
calculations are performed in the same gauge. We will use an on-shell basis for the potential 
where the term in Eq. (^j) does not occur. The matching coefficients for the on-shell potential 
are gauge independent, since scattering amplitudes are measurable quantities. The difference 
between the on-shell matching potential and the Breit Hamiltonian is compensated for by 
a corresponding change in the l/|p' — p| potential, as discussed in more detail in section III 



(see also Refs. |2f^9[). 

In section IV we extend all of our results for the quark-antiquark potential to the quark- 
quark potential and in section V we give the QED limit of our results. Section VI gives our 
conclusions. 



II. MATCHING THE POTENTIAL TO 0{v 2 ) 

The vNRQCD effective Lagrangian has the form 

C = C u + C s + C p . (7) 

The ultrasoft Lagrangian C u involves the fields ip p which annihilate a quark, Xp which anni- 
hilate an antiquark, and A 11 which annihilate and create ultrasoft gluons. The potential La- 
grangian C p contains operators with four or more quark fields including the quark-antiquark 
potential. Finally the soft Lagrangian C s contains all terms that involve soft fields which 
have energy and momenta of order mv. Heavy quarks with soft energy and momenta are 



5 



off-shell and are therefore integrated out, so they do not appear explicitly in C s . The terms 
we need in the ultrasoft Lagrangian are 

c =-l F i» F +v^L jD o_ (p-^ d ) 2 + _p1L 



p 



3m 3 



>Xp- (8) 



The covariant derivative on ip p and Xp contain the color matrices T A and T A for the 3 and 
3 representations, respectively. Here involves only the ultrasoft gluon fields: _D M = + 
ignfjA^ = (D , -D), so D° = d a +igy\jA , D = V-igy^A. The factors of the ultrasoft scale 
parameter \x\j are included to make g = g(fiu) dimensionless in d = 4 — 2e dimensions. The 
ultrasoft gluon field is order v 2 and has dimension 1 - e, so # ~ (mw 2 ) 1_<E . Consistency 
of the v power counting for the covariant derivative in d dimensions then requires nu = rnv 2 
where v ~ v. This reproduces the dependence of \x\j on the subtraction velocity v given in 
Ref. [16|. The terms we need in the soft Lagrangian are 

= E { \?K - 9Xf + ^ M + c q q\ ) (9) 

P,p',9,<2' 

+ V cjyw v P + (V t b z& v P ) ( < /v7 / ^~ B ^>g) j + (V> - x, r - f ) , 

where A g , c q , and y> 9 are soft gluons, ghosts, and massless quarks. The functions U, W, Y, 
and Z are given in Appendix After integrating out the soft quarks the Lagrangian C s 
is no longer manifestly gauge invariant with respect to gauge transformations at the scale 
mv. Therefore, determining the dependence of the soft scale parameter fj,g on v may seem 
more difficult than the ultrasoft case. However, prior to integrating out the soft quarks the 
combination gy s A q is from a covariant derivative, and A q ~ {mv yielding fi$ ~ "mv in 
agreement with Ref. [16|. In Eq. @, g = g(fis)- In general it is important to realize that 
the v scaling of \iy and fig are different. If the matching calculation is performed at the 
scaleQ m, then for this computation \x = ys = — m an d u — 1 (where the usual QCD 
scale parameter is denoted by fi). Therefore, for the matching at m it is not essential to 
distinguish between fig and \x\j. 



Tt is not necessary to match exactly at m. If the matching scale is y — y^ ~ m, then 
one still sets yu = Us = l^h an d v — 1, and factors of \n([i h /m) appear in the matching 
coefficients. For convenience we choose /id = m in this paper. 
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The potential interaction relevant for our calculation is 



£p = ~ V *^r (P, P') ^p'L Vfe/J X-p'[ X-pr- 

p,p' 



(10) 



The factor of /i| e is included so that in d = 4 — 2e dimensions the potential V has dimension 
—2. a, /3, A, r denote color and spin indices. We will use the basis in which the potential is 
written as a linear combination of 1 ® 1 and T a g) T a in color space. One can convert to the 
color singlet and octet potential using the linear transformation 



^singlet 




►octet 





-c f 



1 \Ca — Cp 



n.®i 



(11) 



where C F = (iV 6 2 - 1)/(2JV C ) and C A = N c . We will also need the invariants C x = (iV 2 
l)/(4iV 2 ) and C d = N c - A/N c . These arise in the identities 



T A T B q rpArpB = (C A + C d )T A ®T A + dl®l 



~{C A -C d )T A ®f A + C l l® 1 



T A T B ® T B T A = 
Written as a matrix, the order v ~ 2 Coulomb potential is 



(12) 



y(- 2 ) = (T A <g> f M c 



+ (1 



k 2 V - / k 2 > 

where k = p' — p, and the coefficients VjF^ and have an expansion in a 
The order v° potential includes 

;( T ) M(T) („2 , n /2\ y(T) 



(13) 



= (T A ® T A ) 



+ (1® 1) 



, vf } (p 2 +p /2 , 



2 m 2 k 2 



V (T) V 

Va a(p',p) + H 



(T) 



T(k) 



77T 



where 



A(p',p) 



■ S'(p'xp) 
k 2 



T(k) 



(14) 



3k (Tik (T2 

o-i<r 2 p , (15) 



and (Ti/2 and cr 2 /2 are the spin-operators on the quark and anti-quark. Note that on-shell 
p' 2 = p 2 , but we have written p 2 + p' 2 in Eq. (|TJj) so that the Lagrangian L v is hermitian. 
The tree level diagram in Fig. |l|a generates terms of 0(v 2k a s ), k > —1, in the QCD potential. 
Matching at jj, — m, v — 1 gives 
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i 

a 
a 

a 1 

1 
a 

_3_ 



(a) (b) 
FIG. 1. QCD diagrams for tree level matching. 



vf) 



47ra s (m) 
AnaAm) 



vP = o. 



VP 



vP 



-67ia s (m) 



V, 



(T) 







v. 



(T) 



0. 



(1) 



47ra s (m) , 
7ra<,(m) 
" 3 
= 0. 



(16) 



The annihilation diagram in Fig. |I]b generates terms of order a s v 2k , k > in the potential. 
Using Fierz identities and charge conjugation, these operators can be transformed into the 
basis in Eq. ([TJP and give additional contributions to the matching. Only Vp 1 ^ receive 
non-zero annihilation contributions at tree level: 



V2? = ^-7ra.H, V$ 



(N* - 1) 

2A„ 2 



7r a s (m) 



(17) 



The complete tree level matching is given by adding the terms in Eqs. (16) and (p|) . We have 
found it convenient to distinguish the annihilation contributions by including an additional 
subscript a on the coefficients they generate. The leading-log values of the order v° potentials 
Jl4]) were calculated in Refs. |3yj2^], but are not needed here. Nonzero values for 



in Eq. 

V2 ' are generated in the renormalization group flow below the scale m pOfl , as well as by 
the one-loop matching as we will see below. 

At one-loop the matching onto QCD gives order 1/v terms of the form 

■2 



it 



m Ikl 



f A ) + vP(i 



(18) 



(T 1) 

where the coefficients ' are dimensionless. In a = 4 — 2e dimensions the one-loop 
matching produces a potential with the dependence /i 2<E /|k| 1+2e . We have chosen to define 
V^" 1 -* by taking the d — > 4 limit, which differs from the definition of this operator used in 
Ref . |I . 



A. The QCD Calculation 

To perform the potential matching calculation, we consider the on-shell QQ scattering 
amplitude in QCD and in vNRQCD to order ot 2 s v° . We will use Feynman gauge, regulate 
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infrared divergences with a finite gluon mass^ A, and renormalize ultraviolet divergences with 
dimensional regularization and the MS scheme. Since the calculation is performed on-shell 
the resulting matching coefficients will be gauge independent. 

We begin by considering the QCD diagrams. The most complicated diagram is the QCD 
box diagram [|32| , |33f which includes contributions of order l/v 3 , l/v 2 , l/v and v°, as well as 
higher order terms which we do not need in this paper. 



2o£ ^ fAfB 

fir 



2imn , /k 2 \ /k 2 \ . . 

— '"(aj)- 41 "(a5) (19) 



mp(2p+fc)l 2 4 " v 6 / 12 A(2p + k) 

|A; 2 - 4 p 2 + (5A; 2 - 12p 2 ) In (^) - (4p 2 + k 2 ) In (^) | 



mp (4p 2 — A; 2 ) 
| z7rfc 2 S 2 f,o , , r; , , iL 

3mp (4p 2 — fc 2 ) 
+ ™k 2 T | , . , 2 , 2] _ 



+ 



6mp (4p 2 — A; 2 ) 

2 mP( v-^) ( 4p2fc2 - fc4 ^ 2(4p2+fc2)in (T) 



/ 4 2,2 ,4x, 6A; 2 8A; 2 S 2 k 2 T 

+ (80p 4 - lQp 2 k 2 + k ) ln (-) - — + — — - — 

ln /Ax |56^ _ 12fc 2 A _ 8A; 2 S 2 _ 2k 2 T \ | t / fc \f 4fc2 _ 2A; 2 S 2 
^/c^l3m 2 m 2 3m 2 3m 2 J ^m^lm 2 m 2 

where k — |k|, p = |p| = |p'| and 

R=(p + p')-^i(p + p')-^2. (20) 

The real part of the order l/v amplitude agrees with Ref. [^2] in the limit p — > fe/2. We 
have kept the full p dependence since taking this limit is not justified by the power counting. 
We have also kept imaginary terms generated by the cut amplitude to emphasize how these 
terms are correctly reproduced in the effective theory. The real part of the spin dependent 
order v° amplitude agree with the result in Ref. ||34|| . The real part of the spin independent 
order v° amplitude agrees with Ref. |23|, except for the order v° ln(k) and ln(A) dependence. 



2 Using a finite gluon mass is dangerous in the presence of diagrams with the non-abelian 
gluon vertices. All such diagrams we require here are IR finite in Feynman gauge. 
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The difference is due to the condition |p' + p| <C k which was imposed in Refs. but 
which violates the v power counting. (For the crossed box given below in Eq. (|2~T|a), the order 
v° ln(A) and ln(fc) also differs from Ref. pB[ .) The remaining direct scattering diagrams are 
less complicated since they have no cuts: 



+ 



+ 



+ 



+ perms 



(T A T B ® T B T M 



\. ( k\ n 2 k 



in 



(21a) 



+ 



A\ r 8A; 2 56p 2 12fc 2 A 8£; 2 S 2 2£; 2 T 



+ ]»(*) 

V 77V 



3m 2 3m 2 /// 
6k 2 2k 2 S 2 



+ 



+ 



+ 



- 3m 2 3m 2 

2fc 2 2A; 2 S 2 fc 2 T 
+ 



+ 



' A; 2 



+ ln(*){ 



k^ (4k 2 



31n(^)+4 + 



m z 



+ 



ln^l^ 



4fc 2 S 2 
3m 2 
9A; 2 A 
m 2 



-m 2 
k 2 T 
3m 2 
2k 2 S 2 
m 2 



3m 2 
vr 2 A; 



3m 2 



2m 



4A; 2 A 



k 2 T^ 

2m? ) 



Ap 2 _ 8k 2 A _ 2k 2 S 2 _ k 2 T 

2m? 

f A ) 
2k 2 A 



77T 



(21b) 



--±(2C F -C A )(T J 



21n(4)+ln(4)+4 M ° 



k 2 



2A; 2 S 



2o2 



k 2 T 



3m 2 6m 2 



ia, 
k 2 



2C F (T J 



C A (T' 



f-M-) 

6m z 



(21c) 

(21d) 
(21e) 



and the light quark loops for rif flavors gives: 

6 = g»^«f*){§ln(£) + £}^, (21f) 

while the heavy quark fermion loop gives: 



10 



k 2 



T F {T A ®T J 



M°- 



4k 2 1 
15m 2 J 



In Eq. fl2lD the matrix element 



1 + 



k 2 S 2 3k 2 A k 2 T 



(21g) 



(22) 



3m 2 2m 2 12m 2 

Note that we disagree with Ref. [23[] on the order v° spin independent part of Eq. (|2l|b). 
Ref. |23] has an additional non- logarithmic 1/m 2 term. The difference arises because we 
find a different value for the order k 2 term in the non-Abelian F\ form-factor. Eqs. (pile) 



through (^Tg) agree with Ref. |[23|| . 

The diagrams in Eq. (|l^) and (|2l|) have contributions from several different scales. In 



particular, in the language of the threshold expansion [13], the hard regime gives ln(///m)'s, 
the soft regime gives ln(/i/fc)'s, and the ultrasoft regime gives ln(/i/A)'s. There are also 
iln(A/p) and iln(\/k) terms from the Coulomb divergence in the potential regime. In 
addition to logarithms, all regimes can give constant factors. In the effective theory, terms 
from the hard regime are absorbed into matching coefficients such as the four-quark potential 
operators and the remaining terms correspond to graphs involving modes in the effective 
theory. These graphs are discussed in more detail below. 

Next consider the one-loop annihilation diagrams in QCD. Since the intermediate gluons 
are hard we expect these graphs to include a factor of 1/m 2 , thus giving hard order v° 



contributions to the potential. However, the graph in Eq. 
a factor of m/p, which are order 1/v contributions. 



+ 



to; 
m 2 

:2 



>Cx C d (N 2 - 
^N r 8N 2 



also has terms enhanced by 



(l®l) + (-^-2C 1 )(T^f A ) 



V4iV, 



S 2 -2)(m + 2-21n2) - C A 



(N? 



2N 2 



l) + J_(T A 



x S : 



Oronnnr 



+ 



x S : 



m 
112 

2 

£ 

2 

3 



1 
6 



In 2 



6 



-In (A) 

v m' 



(23a) 



C A \ (T y 



T )—— + (1 



12 W 



+ 



21n2 



W - 1) 
2N r 



§. ymnswir 



+ 



Am 2 



(2C p — C/ 



3 

(T A ®T 







y 


}• 




1 







(23b) 



+ (1 



(N 2 



2N r 



x S 2 + hi(-H -4 + 21n(-^-) +41n( — ) , 23c 

L p p v A ' x m' v m y J 



11 




+ perms 



C F {(T A ®f A )^- + (l®l) 



A 



S 2 hn (^) + In (£) + 2 



xS 2 (ln(4)-21n2 + ^ + Z7rl 
L v m z/ 15 J 

and the light quark loop for n/ flavors gives: 



1 

'iV c 
31 



(i Vc 2 ~ 1) 

2N r . 



2N r . 



(23d) 



(23e) 



ia 2 s n f T F ( A 



while the heavy quark loop gives: 



i n ix (iv c 2 -i; 



(T <g> T ) — — h (1 ® 1) 



x In ^— -21n2 + - + Z7T 



2AT r 



(23f) 



^ ^ (T* 



x S 2 



f A )-J- + (l®l^ iVc2 



AT, 



In 



n? 



■i) + 3 



2A, 



(23f 



In the limit p — > fc/2, the results in Eq. (|23|) agree with Ref. P2| , except for Fig. (P3|c) 
and differences that can be accounted for by the fact that we are using MS rather than 

~~|c), Ref. p2] has a ir 2 /k term which for us is 



an on-shell subtraction scheme. For Fig. 
— 7r 2 /(2p). This sign for the tt 2 /p term is in agreement with Ref. ||35|| . The imaginary parts 
of these one-loop annihilation amplitudes also agree with Ref. ||. Note that in Eqs. 
(0), and (H), a s = a s {n). 



B. The Effective Theory Calculation 

The effective theory contains potential, soft and ultrasoft loops. We have organized the 
terms by their order in the velocity expansion. The order in v can be determined using the 
v power counting formulaf] in Eq. (40) of Ref. [IB]. A loop graph with two insertions of the 
Coulomb potential contributes to a 1/v 3 potential, and with an insertion of one Coulomb 
and one potential contributes to a l/v potential. A soft loop with two vertices of order 



3 Here the power of v is given for the amputated diagram, so unlike Eq. (40) in Ref. [IB 
the factors of v associated with external lines are not included. 
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a and a' contributes to the v a+CT ' 2 potential. Graphs involving the exchange of an ultrasoft 
gluon begin to contribute to the potential at order v°, etc. 



1. Order 1/f 3 

In the effective theory taking two insertions of the tree level Coulomb potential in a loop 
gives the only diagram that is order a 2 /v 3 : 

XX -C^™^!?^)- m 

Taking /i = m and v — 1 and using the tree level value of Vj T \ this graph exactly reproduces 
the order a 2 /v 3 "Coulomb singularity" term in the QCD box diagram in Eq. (fOJ). Thus, 
there is no matching correction at this order. 



2. Order l/v 2 



At order a 2 /v 2 in the effective theory, the only non-zero diagram involves the exchange 
of soft gluons, ghosts and quarks: 




k 2 



-X 



4T F n f - HCy 



{Xd)} 



+ 



20T F n f - 31C A 
9 



(25) 



where 1/e = 1/e — 7 + ln(47r) is the combination subtracted in MS. The divergence in this 
graph is responsible for the one- loop running of VX . At one- loop, to order v the effective 
theory counterterms (and running of the potential) were computed in Ref. pDJ and from 
now on the 1/e dependence will be dropped. The sum of order a 2 /v 2 terms from the QCD 
diagrams in Eqs. (|T9|) and (^) is 

X(X 



k 2 



-X 



AT, 



4T F n f -UC A /X 20T F n f -31C A 
3 n VX + 9 



(26) 



At the scale fi = us = m the MS values of Eq. (g^) and Eq. (p6|) are identical, so there 
is no one- loop matching correction to V c . The difference in ln(/i)'s in Eqs. ( p6|) and ( p5|) is 
the change in /3-function from rif + 1 to rif flavors. We would expect a new contribution 
to V c only if the QCD graphs have a contribution from the hard regime where energy and 
momenta are order m. In Eqs. (|2ljb), (|2l|c), and (|2l|d), the factors of —3 ln(m 2 ) +4 are from 
the hard regime, but they sum to zero. For the ln(m) terms this is a consequence of the 
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Ward identity derived in Ref. ||. The constant factor (20Tjrn/ — 31Ca)/9 vanishes in the 
matching condition at m, and is carried to scales below m by the soft modes. 

Note that if at the scale v = vj. ~ k/m we integrate out the soft modes, then the 
Coulomb potential for the theory below this scale, V c , will obtain an additional contribution 
from this second stage of matching: V c (vk) = Vc{ v k) ~ (20Tpn/ — 31CA)oi 2 (mvk) /9. In 
this expression V^ T '(vk) is the value of obtained from running this coefficient from v = 1 
to v — Vk using its two-loop anomalous dimension. This reproduces the constant factor that 
is typically associated with the Coulomb potential at next-to-leading order (see for instance, 
Ref. [|H]). Similarly, one can obtain the additional terms from the matching contributions 
for the v° potentials at mv from the value of the soft loops given in Eq. ( pi]) . 



3. Order 1 jv 

The possible order ot 2 s /v diagrams with soft gluons vanish explicitly, so the only order 
a 2 /v diagrams in the effective theory involve two iterations of the potential. There are two 
diagrams with insertions of V^ 7 " 1 : 

V c v- V c V c V c •f\;(T)l2 

X^jK^+ AE^^y^ = i -^±T A T B ®f A f B (l F + 2p 2 Jo), (27) 

where the integrals Iq and Ip are given in Appendix 0. The dependence on V^ 1 -* is not 
needed since tree level matching gives V^ 1 -* = 0. In the first diagram in Eq. fl2"7|), the cross 
denotes an insertion of the p 4 /m 3 operator from Eq. @. The second diagram is nominally 
order a 2 s /v 3 ; however it depends on the heavy quark energy E = p 2 / (2m) — p 4 / (8m 3 ) + . . .. 
When the energy is expanded in terms of momenta, the graph includes a contribution of 
order a 2 /v which we indicate by the pre-factor AE in Eq. fl27D . Each of the diagrams in 
Eq. (p7|) has an IR divergence that is not regulated by A, but the IR divergence in the sum 
of the integrands for the two diagrams is regulated. 

Additional a 2 /v diagrams are generated by including one insertion of the Coulomb po- 
tential and one order v° potential: 



V c _ V V _ V c 

+ 



iVPVP rpArpB ^ rpArpB ( J p + 2 J? I Q ) 

m v ' 



(28a) 



v c v s v s v c 



+ 



Vc Vs T A T B ® f A f B (21 p S 2 ) 



(28b) 
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V A Va _ V c ;v(T)v( T ) 

+ = m TATB ® ( 2P A(P '' P) /a ) ' (28C) 



+ V^V" = tkc ^ T A T B ® T A T B 



2m 

4<n • a 2 (I P - 3I B ) + R(12l£, - 12I A + 3/ 



k-a 1 k-a 2 (12/ c + 3/ ) 
(28d) 




iV (T) V (T) , 

' T A T B ® T A T B (S 2 2/ P ) 



m 



rriA 



T <g> T ( S 2Jp) . (28e) 



m 

The last two diagrams involve insertions of terms in the potential generated by the tree 



level annihilation diagram. In Eqs. (|2T|) and (|28|) the dependence on 1 <g> 1 potentials whose 



coefficients vanish at tree level are not shown. Thus, both the 1 ® 1 and T ' ®T contributions 
are only shown in Eq. (p8|e). The new integrals Ip, I a, Ib, Ici an d Id that appear in Eq. ( |28| ) 
are given in Appendix |B|, and R is defined in Eq. ( PD|) . 

The infrared divergences in the 1/v amplitude are due to the Coulomb singularity. The 
divergences in the full theory are reproduced by the potential loops in Eqs. ([27]) and (]28|). 
The imaginary terms in the amplitude from the QCD box graph exactly agree with those 
in the effective theory, and do not contribute in the matching coefficients, as expected. 
Subtracting the effective theory graphs in Eqs. (^7]) and (|28|) from the order ot 2 s /v terms 
in Eq. ([19|) , (pl]) a,b and fl23"l) c gives the matching result for the V^ -1 -* potential at [i = m, 
u=l: 

VP = cUm){^-%), Xf=oJ(m)^. (29) 

For the color singlet channel this gives the matching coefficient = a 2 s (m){C F /2 — CfCa) 
in agreement with Ref. f23|. Note that the real part of the 1/v amplitudes in the full and 



effective theory have a complicated dependence on p and k, but the momentum dependence 
of the matching for the direct potentials in Eq. (p9|) is only of the form 1/ |k|. For the 
annihilation graphs the 1/p terms cancel between the full and effective theories. There 
would have been a 1/p matching coefficient for the annihilation graphs if the results of 



Ref. [f22j had been used for the full theory graph. 

The matching coefficients in Eq. (|29| ) correspond to a contribution to the 1/v potential at 
the scale fi = m, which does not come from the hard part of any graph. This is in apparent 
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contradiction with the threshold expansion. It also appears to disagree with Refs. [19] and 



53J , where the l/|k| four-quark potential operator is said to only arise at the scale mv. 
However, in Refs. [19| and a formulation of NRQCD is being used in which off-shell 
potential field components have not been integrated out for scales mv < fi < m, but instead 
are considered to be dynamical fields in the Lagrangain (see Refs. [0,|1O[). In the full theory 
the l/|k| terms come from three types of graphs, shown in Eq. flT9| ) and Eqs. (]2"T|a,b). At 
the scale m there are now effective theory graphs analogous to those in Eq. (|2l]b) but with 
two A potential gluons and one A 1 potential gluon, which reproduce the l/|k| term. For 
the box and crossed box in Feynman gauge the l/|k| terms are reproduced by contributions 
that can be associated with the potential momentum regime. Thus, in Refs. and [31 



the l/|k| potential in Eq. ( p9|) effectively exists at the scale m. In our approach off-shell 
potential gluons and soft quarks are integrated out, so the matching for the l/|k| potential 
is not simply given by the hard part of the QCD diagrams. 



4- Order v c 



For the order a 2 v matching we will consider the direct and annihilation diagrams sepa- 
rately. The sum of the order a 2 v° terms in the direct scattering QCD diagrams in Eqs. ([19]) 
and (0) is 

(l®l)cj4-2S 2 + 21n( * 



rn 

,2 



(1 ® l)C x 
(T A ® f A 



p 2 j31C A 1QC A 
k 2 \ 9 3 

2Cf — 3C A — Cd 
7C d A3C A 







5 In 




3 


(i). 




3 


(-) 


+ - 


v 





\Cf 2Cd 



k 

2C A ) In 



+s 2 

+T< 



6 6 
31 C A 



. /8C F 31C 



6 



6 



4CW^ln@-4^1n(9} 



Cd 
2 



17 C A 4CV C A 



54 



9 W 



7C* 



49^4 
108 



Cf _ 5Ca , //A _ Ca j n 1 
3 18~ 11 W T 11 Wj 



11 W/\fc2 ~ "3 2 12 



p 2 S 2 _ 3A T 
^2 _ 1 2 ~ 12 
S 2 3A T 1 4T; " 
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(30) 
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The effective theory contribution from order a 2 v diagrams with soft gluons, ghosts or quarks 
is: 



1) 



14d 



In - 



+ 



s i 



mf 
7C d 



3 

43 



(31) 



38p 2 



7A H 1 

9 18 



6 ~ ^ k > 12 ' 3 



+n f T F 



ln( 



31CV 7C A A 



9fc 2 



6 



14C A S 2 13C A T 



27 



108 



A; 2 ^ + 9 



,2 

A- 2 



12 



There are also non-zero order a 2 t> diagrams with an ultrasoft gluon and one insertion of 
the Coulomb potential. These diagrams were calculated in Ref. [pOj 




3m 2 



A > 



Ci(l ® 1) 



a s {ns)a s {nu) In (— ^ 



A >2 



(32) 



In Eq. ( [32]) the ellipses denote all possible diagrams that are generated by attaching the 
ultrasoft gluon to two fermion legs. Subtracting the sum of Eq. (|3T|) and (|32|) from Eq. ( |30| ) 
and setting /i = [is = flu = m gives the one-loop matching contribution for the order v° 
direct potentials. At order v we find that all the infrared divergences in QCD are matched 
by infrared divergences from the effective theory graphs with ultrasoft gluons. All ln(fc)'s 
in the full theory are matched by ln(fc)'s in Eq. fl3~T|). The contributions to the matching 
coefficients at one-loop are: 







vf) = , V« = 

VP = -A(C F + C A )a 2 s (m), 
VP = (\C d - \C A - \C F ) a 2 (m) , 
VP = -l(C F + C A )a 2 (m) , 

VP = (2C F - %C d - \C A + ±T F ) a 2 (m) , V 



y« = -2d a 2 s (m) 



(33) 



VP 



(i) 



, 

f G x a 2 (m) 



The total order a 2 v matching coefficients are given by adding the tree level values in Eq. ([16]) 
to the results in Eq. (|3"3"D, and are summarized in Table | at the end of the paper. 

Next consider the matching of the order a 2 v annihilation contributions. Since there are 
no corresponding diagrams in the effective theory, the sum of the a 2 v terms in Eq. ( |23|) 
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directly give the matching coefficients. This sum is infrared finite. Matching at // = m, 
v — 1 we find that the one-loop annihilation contributions to the potential coefficients are: 

Vg? = {%- 2Ci) (ztt + 2 - 2 In 2) a 2 s (m) 



v 4^ 



/ 7T i — > a 2 Jm) 



iV c l 36 F 3 V 3 

d , d(iV 2 -l) 



s 1 



V« = {f + V + 2 - 2 ln2) a2(m) 



(iV 2 - 1) f 109^ ^ n f T F / , 5 x 87* 



2iV c 2 I 36 3 
vff = -2(^ - 2d) (ztt + 2 - 2 In 2) a 2 (m 



4iV c 

Vg = -2{| + + 2 - 21n2) a2(m) . (34) 

The imaginary terms in these potentials contribute to the cross section for annihilation of 
a color octet heavy quark and anti-quark into light hadrons, and agree with the results of 
Ref . || . The total annihilation contribution to the order a 2 v° matching coefficients are given 
by adding the tree level results in Eq. ( [T7D to the results in Eq. ([54]) , and are summarized 
in Table |. 

If a different matching scale, fit, had been used then the coefficients in Eqs. (|33| ) and 
(51) would also depend on ln(/^/j/m). Since the prediction for observables is independent of 
Hh the most convenient choice, fit = m, has been adopted. 

In the threshold expansion, the full QCD diagram is the sum of hard, soft, ultrasoft and 
potential graphs. The soft, ultrasoft and potential graphs are the graphs in the effective 
theory up to renormalization effects such as the dependence on the scales us and The 
matching conditions in Eq. (|33|) and (0) can also be computed from the hard part of the 
full theory graphs, and we have verified that this gives the same result as in Eqs. (|33D and 



(34). The matching onto four-quark operators was first considered by Pineda and Soto in 
the context of pNRQCD [|12|| . In their approach, the hard parts of the vertex correction 
and wavefunction renormalization are included in matching coefficients in the single heavy 
quark sector of the Lagrangian. The direct matching for four-quark operators is given by 
the hard part of the box and crossed-box. The hard part of our box and crossed-box agree 
with Ref. |[I2| , except for the finite part of the o\ ■ a 2 terms, which causes our and 
VP to disa gree with theirs. This is related to the treatment of epsilon tensors and the 
non-relativistic reduction of matrix elements of spin operators in d dimensions. We have 
chosen to take the lowest order term in the matrix element of ^ a >y a ~'y 13 } <gi 7[ ce 7 T 7 i g] to be 
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e a @ ak e a /3 T k' <7\C?2 ) an d have used the 't Hooft-Veltmann scheme for the epsilon tensors so 
that and e y ' fc are only non-zero when the indices are in four and three dimensions 

respectively. This convention was used in both the soft loop integrals as well as when cross 
checking our result by computing the hard part of the box diagrams using the threshold 



expansion. This scheme dependence is related to the issue of evanescent operators |36- 38 



The annihilation results in Eq. ( p4|) agree completely with Ref. [[12 



Note that in a leading-log expansion the two-loop anomalous dimension is needed at the 
same time as the one- loop matching results. In the color-singlet channel the 1/v 2 two-loop 
anomalous dimensions is known |l],@], since the running of the Coulomb potential is still given 
by the QCD /^-function at this order. The result for the two-loop anomalous dimension for 
V^ 1 ' will be presented in another publication [R3]. 



III. TERMS IN THE POTENTIAL WHICH VANISH ON-SHELL. 



For the scattering Q{p\, p) + Qip®, — p) — ► QiPzi P') + Q(p% — P') consider adding 



y(0) 



Vl? (T A ®T A )+V A ^(1®1 



+ 



vP 



12 



p 



2\2 



4m 2 k 4 



(35) 



to the potential in Eq. (P~4[). Here k = p' — p and by energy conservation p® — p\ = p\ — p\. 
On-shell the potentials in Eq. fl3"5|) vanish, since p\ = p®, p® = p®, and p 2 = p' 2 . However, if 
we work off-shell, they are valid terms and in fact show up in many calculations that make 
use of time-ordered perturbation theory. 

Matching to the tree level diagrams in Fig. [l] with p\ 7^ p\ gives a contribution to the 
potentials in Eq. ([35|). In Feynman gauge one finds 



y(T) 



K A2 







4vr« s (m), V$ = 0, 
while in Coulomb gauge one gets a different answer 

0, V$ = 0, VS = -4na s (m 



V A\ 



V (l) 
V A2 



V (l) 
K A2 







0. 



(36) 



(37) 



Unlike the on-shell potentials discussed in the previous section, the matching conditions 
for off-shell potentials are gauge dependent. Using the expressions for the soft vertices in 
Feynman gauge in Appendix [A], the one loop anomalous dimensions are 

„ d V (T) _ nn ( s2 d V (T) _ 9 (l) d (i) 



0. 



(38) 
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Including the potential in Eq. ([35]) modifies the matching condition for the l/|k| poten- 
tials and also makes this matching gauge dependent. This is because in the effective theory 
there are now two new order ot s /v diagrams: 

(39) 

In Eq. (^) the box labeled by Va denotes an insertion of both operators in Eq. fl3"5"|), the 
ellipses denote operators that vanish by the equations of motion (proportional to p\ — p 2 /2m, 
P3 ~ p /2 /2m etc.), and we have defined 

V1 T) = Vff + VS , V%> = Vffl + VS ■ (40) 
In Feynman gauge the matching conditions in Eq. (]29D now become 

V£? = a2(m)(^-^), V«=a2(m)C ll (41) 

while in Coulomb gauge 

VI T J = C A a 2 s (m), Vjg = 0. (42) 

The l/|k| potential is often referred to as the non-abelian potential. From Eqs. fl4"I|) and 
(0), we see that the l/|k| potential only vanishes in QED if the potential is taken to include 
off-shell components and Coulomb gauge is used. 

The result in Eq. (]3Qf) shows that in the off-shell matching potential one can make the 
replacements 

VP - VP + , (43) 
VP - + C (1) , 



vP - vP + -^y c (T) 

fe ft 327r 2 c 



C (1) -j(^ + Q)C (T) 



- vP + ^- 2 vp c lC <*> , 

for arbitrary C^' 1 ^. For instance, at the matching scale taking = 0, P^ = — 8ira s (m) 
effectively transforms the Feynman gauge result in Eqs. (|36| ) and fl4T| ) into the Coulomb gauge 
result in Eqs. ( |37| ) and fl42|). Furthermore, taking = and C ( - T ' ) = 4:%a s (m) transforms the 
off-shell Coulomb gauge result into the on-shell result in section II. These transformations 
convert terms in the potential that are order a s v to order a 2 s /v. Similar transformations for 



the position space color singlet potentials have been pointed out previously in Refs. 25-27 
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IV. THE QUARK-QUARK POTENTIAL 



The quark-quark potentials in the effective theory can be defined in the same way as the 
quark- ant iquark potentials. The only difference is that the terms are now the coefficient 
of the T A <g> T A tensor, rather than the T A ® T A tensor. The computation of the quark- 
quark potential is almost identical to that of the quark-antiquark potential. The result can 
be obtained from the quark-antiquark potential by omitting the annihilation terms, and 
making the substitutions Cd — > —Cd and T A — > T A . The change in sign of Cd arises from 
the identities in Eq. JI2]) . Replacing T by T in these equations changes the sign of the Cd 
terms. 

The potential in the symmetric and antisymmetric QQ color channels (the 6 and 3 for 
SU{2>)) are given by 



^symmetric 




1 


iVc-l 

2N C 




Vim 


^antisymmetric 




1 


N c +1 
2N C 




Vt®t 



(44) 



The spin-1 QQ combination is spin-symmetric and the spin-0 QQ combination is spin- 
antisymmetric. For identical fermions in the initial state, for the symmetric spin-color states 
we must antisymmetrize the potential in the momenta, V = V(p, p') — V(— p, p'), and for 
the antisymmetric spin-color states we must symmetrize the potential in the momenta, 
V = V(p, p') + V(— p, p'). This corresponds to including the crossed diagrams in the 
computation of the potentials. 



V. QED 

It is straightforward to obtain the QED potential from our results. For oppositely charged 
particles of charge ±Q, the QED direct potential is given by Q 2 (^i®i — ^r®?)> where V^i 
and V T(S) f are given by our QCD results with C\ = Cp — T F — 1, Ca = Cd = 0, and a s — > a. 
The on-shell l/|k| potential does not vanish in this limit. In the results for the annihilation 
potentials at the scale m in Eqs. (|T7|) and (0), explicit factors of A^ c were included in the 
Fierz transformation, so it is simplest to just separately list the QED limit of these results: 

yCn = _(f7r + 2 -2 In 2) a 2 (m)Q\ 

Vg = vr a(m) + { - f + ^(2 In 2 - | - z^r) } a\m)Q 2 , 
V£ } = 2(ztt + 2 - 2 In 2) a 2 (m)Q 2 , 

= . (45) 
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For e + e~ we have n/ = and the terms in our v° potentials agree with Ref. [JT5 

For identical particles with charge Q, there is only the direct potential contribution, 
which is given by Q 2 (Vi®i + Vt®t), with G\ = Cp = T F = 1, Ca = = 0, and a s — > a 
as above. As discussed in Section [TV], including the crossed diagrams gives a final potential 
that is symmetric or antisymmetric in the momenta depending on the symmetry of the spin 
state. 



VI. CONCLUSION 

We have computed the QQ and QQ scattering amplitudes in QCD to order v 2 , and 
compared our results with previous calculations. We have also computed the scattering 
graphs in vNRQCD, and computed the matching condition at /i = m between the two 
theories. The matching potential was computed using on-shell matching, omitting terms 
which vanish by the equations of motion. The result was compared with approaches that 
include terms in the potential that vanish on-shell. The one-loop matching coefficients are 
summarized in Table |, and can be combined with the two-loop running to give the potential 
at next to leading log order. The computation of the heavy quark production current at 
next to leading logarithmic order uses these results and will be discussed in a subsequent 
publication 
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TABLE I. Matching coefficients for the quark-antiquark potential at fi = m, v = 1. The 
tree-level contributions are the order a s terms, and the one-loop corrections are the order a 2 
terms. The values are for the on-shell potential, so all off-shell potentials such as the Va potentials 



in Eq. (35) are set to zero. Contributions to the matching from annihilation diagrams are given 



separately and are denoted by an extra subscript a. 



v (l) 

v 2,a 



VP 4vra s (m) 

VP 

Vf al{m)(lC A -\C d 

VP Aixa s (m) 

V« 

VP -Q7ra s (m) - A(C F + C A )a 2 (m) 

VP 

VP -§vra s (m) + (\C d - \C A - ±C F )a 2 s {m) 

V« -2d a 2 s {m) 

VP -§7ra s (m) - \{C F + C A )a 2 (m) 

VP 

VP (2C F - %C d - \C A + ±T F )a 2 (m) 

V? ] f C x a 2 (m) 

Vg? + (^-C d - 2d)(i7r + 2 - 21n2) a 2 (m) 

+ i{f -4C J? + 2^(21n2-|-t7r)- ^} a 2 (m) 

V$ ^ 7ra s (m) + + ^C d }(m + 2 - 21n2) a 2 (m) 

+ ^(lf CU - *Cf + ^(21n2 - | - ztt) - ^ j a 2 (m) 



VfJ -2 (^d - 2d) (ztt + 2 - 2 In 2) a 2 



a; m 



- 2 {wpi + ^mf C ^ + 2 - 21n2) a 2 (m) 



23 



APPENDIX A: COEFFICIENTS FOR THE SOFT LAGRANGIAN 



The coefficient functions for the soft Lagrangian in Eq. (|9]) in Feynman gauge are: 
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Am 2 q [ 



[p + p') ic F e ijk cr k (p + p') q (1 - c|)V(p - p' + q) J ' 



Am?q° 



Am 2 q° 
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4 q 2 S ij i5 ij q°c s ia- ■ (p' x p) 2q°6 ij (p' 2 - p 



2\2 



= 



w$ = — 



4m 2 q° 2m 2 (p' — p) 2 4m 2 (p' — p) 4 

(P " P' + q/ 



IP-p 



oo 
y(o) 

y(2) 
z o 

7(2) 



1 

2m 2m(g ' 12 
o 



(i) 



o; 



W. 



(i) 



;o 



2mg° 

-<•/" = q ■ (P + PO + ic F <t ■ [q x (p - p')] 

2m(p' - p) 2 



-q' 

2mg° ' 



h* > = s - 

J 2m 



(P' " P) 2 ' 



Cp q° _ c s icr ■ (p' x p)g° 
8m 2 4m 2 (p' — p) 2 



(P' " P) 2 ' 



(0) 







4m 2 



CD 

2 



C S MT • (p Xp) 



(P' - P) 2 



0, z 



,(2) 



(1) 



(p + p')'-^,f[(p-p') x <x] 
2m(p' - p) 2 



0. 



These expressions differ from those in Ref. [^0| by terms proportional to p' 2 — p 2 , and the 
values in Eq. ( |Al|) are the complete on-shell expressions. The p' 2 — p 2 terms were not 
needed in calculating the one-loop running of the on-shell order v potentials in Ref. p0| . 



In section III, these terms were used to compute the running of the off-shell potential in 



Eq. (38), and this result depends on the fact that we used the on-shell soft Lagrangian. The 
coefficients in Eq. can be written in a manifestly Hermitian form, however we have 

instead used q' = q + p — p' and q' = q° + p° — p' to eliminate q ' since this form is more 



convenient for calculations. Reparameterization invariance [40] can be used to eliminate cs 
by the relation cs = 2cp — 1. The running of Co and cp are given in Refs. 



APPENDIX B: INTEGRALS 



The effective theory integrals that appear in Eqs. (p7|) and (|28f) are 



and 



d 3 q 



d 3 q- 



[(q-p) 2 + A 2 ][(q-p') 2 + A 2 ][q 2 -p 2 -ze] 87r|p|k 2 Vk 2 7 ' 
1 1 



q-p) 2 (q-p') 2 8|k| ' 



d 3 q 



[(q-p') 2 + A 2 ][q 2 -p 2 -ie] 16|p| 8tt| P 



' ln^l p l 



d 3 q 



q' 



[(q " P) 2 + A 2 ][(q - P') 2 + A 2 ][q 2 - p 2 - ie] 



(Bl) 



(B2) 
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d 3 q 



in! 



q'q 



[(q - P) 2 + A 2 ][(q - p') 2 + A 2 ][q 2 - p 2 - ie] 



^I B + (p'-pY(p'-pyi c 



+(p'+pr(p'+ P y/ D , 



where k = p' — p and 

(2|p||k| - k 2 )vr - 2*k 2 In (^) -4*p 2 ln(g) 



Ia = 



In = 



Id = 



167r|p|k 2 (4p 2 - k 2 ) 

(2|p|-|k|)7r + 2»|p|ln(^) 
16tt(4 P 2 - k 2 ) 

(2|p| |k| - k 2 )vr - 2i(4p 2 - k 2 ) - 2ik 2 In _ 4ip 2 In (g) 



1 



327r|p|k 2 (4p 2 - k 2 )^ 



327r|p|k 2 (4p 2 - k 2 ) 
(2|p| - |k|) 2 (4|p| + |k|)|k|7r + 2*k 2 (4p 2 - k 2 ) 



-2*k 2 (12p 2 - k 2 ) In (HM) _ 4 *p 2 (4p 2 + k 2 ) In (^) 



(B3) 
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